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ABSTRACT 

The relation between the 2 — 10 keV, long term, excess variance and AGN black hole 
mass is considered in this work. A significant anti-correlation is found between these 
two quantities in the sense that the excess variance decreases with increasing black 
hole mass. This anti-correlation is consistent with the hypothesis that the 2 — 10 keV 
power spectrum in AGN follows a power law of slope —2 at high frequencies. It then 
flattens to a slope of —1 below a break frequency, Vhfb, until a second break frequency, 
vifb, below which it flattens to a slope of zero. The ratio Vhfb/^ifb is equal to 10 — 30, 
similar to the ratio of the respective frequencies in Cyg X-l. The power spectrum 
amplitude in the (frequency x power) space does not depend on black hole mass. 
Instead it is roughly equal to 0.02 in all objects. The high frequency break decreases 
with increasing black hole mass according to the relation Vhfb= 1-5 x 10~ 6 /(M/10 7 
Mq) Hz, in the case of "classical" Seyfert 1 galaxies. The excess variance of NGC 4051, 
a Narrow Line Seyfert 1 object, is larger than what is expected for its black hole mass 
and X-ray luminosity. This can be explained if its Vhfb is 20 times larger than the 
value expected in the case of a "classical" Seyfert 1 with the same black hole mass. 
Finally, the excess variance vs X-ray luminosity correlation is a byproduct of the 
excess variance vs black hole mass correlation, with AGN accreting at ~ 0.1 — 0.15 
the Eddington limit. These results are consistent with recent results from the power 
spectral analysis of AGN. However, as they are based on data from a few objects only, 
further investigation is necessary to confirm that there is indeed a "universal" power 
spectrum shape in AGN (in the sense that the value of the power spectrum parameters 
of most AGN will be distributed around the "canonical" slope, and amplitude values 
listed above). One way to achieve this is to determine the excess variance vs black hole 
relation more accurately, using data from many more objects. This will be possible in 
the near future, since it is easier to measure the excess variance of archival light curves 
than to estimate their power spectrum. The excess variance vs black hole relation can 
therefore play an important role in the study of the X-ray variability scaling with 
black hole mass in AGN. 
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1 INTRODUCTION 

Since the beginning of the active galactic nuclei (AGN) X- 
ray variability studies it was noticed that more luminous 
sources show "slower" variations. Barr & Mushotzky (1986) 
were the first to show that the "two-folding" time-scale (i.e. 
the time-scale for the emitted flux to change by a factor of 
two) was faster in lower luminosity objects. Later, Lawrence 
& Papadakis (1993) and Green et al. (1993), using the results 
from the power spectral density function (PSD) analysis of 
the EXOSAT "long looks", showed that the PSD ampli- 
tude at a given frequency decreases with increasing source 
luminosity. Nandra et al. (1997) and George et al. (2000), 



using the so called "excess variance" (of ms , i.e. the variance 
of a light curve normalised by its mean squared after cor- 
recting for the experimental noise) found an anti-correlation 
between excess variance and source luminosity. Similar re- 
sults were also presented by Leighly (1999) and Turner et al. 
(1999), who also found that a particular group of AGN, the 
so called "Narrow Line Seyfert 1" (NLS1) galaxies, show sys- 
tematically larger variability amplitude when compared to 
the same luminosity "classical" AGN (i.e. Seyferts with pre- 
dominantly broad permitted lines, BLSls). Further progress 
in the study of the longer-term X-ray variability has been 
afforded by RXTE, which has allowed systematic observa- 



2 I.E. Papadakis 



tions over time-scales longer than before (i.e. months, years). 
Markowitz & Edelson (2001), using 300 day long RXTE light 
curves of nine Seyfert galaxies, estimated their excess vari- 
ance and found a significant anti-correlation between al ms 
and source luminosity even on these long time-scales. 

Most of the previous studies have examined the depen- 
dence of the "variability amplitude" on the source luminos- 
ity. However, during the last few years, the mass of the cen- 
tral black hole has been estimated for many AGN (see e.g. 
Woo & Urry 2002 for a recent compilation of AGN black 
hole mass estimates). It is therefore possible to investigate 
the dependence of the variability amplitude on the black 
hole mass (Mbh), a fundamental property of AGN. Recent 
studies that have followed this approach, have shown that 
there exists a significant anti-correlation between alms and 
Mbh, using ASCA light curves of BLSls and NLSls (Lu & 
Yu, 2001; Bian & Zao, 2003). 

The primary goal of the present work is to investi- 
gate the relation between excess variance and Mbh using 
the al ms data of Markowitz & Edelson (2001), together 
with the recent al ms measurement of a radio-quiet quasar, 
namely PG 0804+761 (Papadakis, Reig, & Nandra, 2003). 
Markowitz & Edelson (2001) have already shown that ol m s 
is strongly anti-correlated with the source luminosity. It is 
important though to see if there exists a correlation between 
arms and Mbh as well. In this way, one can relate the ob- 
served X-ray variations to the physical system itself, since 
Mbh is an important physical parameter of AGN. Based on 
the fact that the integral of the PSD is equal to the variance 
of a light curve, the o-J: ms -MBH relation is used in order to 
examine the intrinsic shape of the PSD in AGN, and how 
this correlates with Mbh- Obviously, a single alms value can 
correspond to different power spectral forms. As a result, the 
^Ls-Mbh relation cannot provide us with a unique answer 
as to what is the intrinsic shape of PSD in AGN, and how 
this might correlate/change with Mbh- However, assuming 
a certain PSD shape, the study of the ov ms -MBH relation 
can be a powerful tool in assessing the reality of the assumed 
PSD shape. 

In this work, the crJ ms -MBH relation is used in order to 
examine the recent results from the power spectrum analy- 
sis of RXTE, XMM-Newton and ASCA light curves of AGN 
(e.g. Uttley, McHardy,& Papadakis, 2002; Markowitz et al., 
2003). These results have increased significantly our knowl- 
edge of the intrinsic PSD shape of AGN. However, high qual- 
ity PSD estimation is currently possible for a few objects 
only, due to the lack of a large number of data sets that 
meet the necessary requirements (i.e. high signal-to-noise, 
long, well sampled light curves). The results presented in 
this work are also based on the use of a small number of 
objects, which have high quality, long term alms measure- 
ments. However, although the number of objects with ac- 
curate PSD estimation may not increase significantly in the 
near future, the study of the of ms -MBH relation can be ex- 
tended soon to incorporate the use of cr rms measurements 
from a large number of objects. For example, the are nu- 
merous AGN observations with RXTE and ASCA, which 
last typically for ~ 1 day (assuming a ~ 30 — 40 ksec ex- 
posure time). Although the PSD estimation is difficult with 
these light curves, the measurement of their a rm s is easier. 
Taking also into account the wealth of the new data that 
are currently provided by XMM-Newton, and CHANDRA, 



it will be possible to define the ffrms"MBH relation more 
accurately in the near future. Consequently, the (tLs'Mbh 
relation will probably play an important role in the study of 
the X-ray variability scaling with Mbh in AGN. 



2 THE DEPENDENCE OF LONG TERM 

VARIABILITY AMPLITUDE ON BH MASS 

Markowitz & Edelson (2001) have calculated the excess vari- 
ance for nine Seyfert galaxies using 300 day long, RXTE 
(PCA), 2 — 10 keV light curves with a uniform 5 day sam- 
pling. Recently, Papadakis et al. (2003) calculated the excess 
variance of PG 0804+761 using a year long, RXTE (PCA), 
2 — 10 keV light curve with a 3 day sampling. The sam- 
pling and duration of the quasar light curve is similar to 
the respective properties of the Seyfert galaxy light curves 
of Markowitz & Edelson. Therefore the quasar a rms mea- 
surement can be used together with the Seyfert values. The 
a^ms values, together with Mbh estimates, are listed in Ta- 
ble 1. In most cases, the Mbh estimates are taken from 
Kaspi et al. (2000), and they correspond to the mean of 
their "rms FWHM" and "mean FWHM" estimates. All ob- 
jects with a a rm s estimate are BLSls except from NGC 4051 
and MCG 6-30-15, which are classified as NLSls. 

In order to compute the uncertainty of the al ms esti- 
mates, the uncertainty of the PG 0804+761 al ms measure- 
ment was estimated first, based on the PSD results of Pa- 
padakis et al. (2003). As al ms is roughly equal to the PSD 
integral from the highest to the lowest sampled frequency, 
using the best-fitting parameter values from the power-law 
model fitting to the PSD, it is straightforward to com- 
pute the uncertainty of the PSD integral over the sampled 
frequencies, and hence alms itself. The result shows that 
[error(alms) / Crms]^ 0.2. The uncertainty of the PSD inte- 
gral depends mainly on the frequency range and frequency 
resolution of the PSD. These properties are determined by 
the number of points and the length of the light curve. Since 
all the alms values listed in Table 1 are estimated by light 
curves with almost identical length and number of points, 
the "signal-to-noise" ratio of the al ms values should be ~ 0.2 
in all cases. The al m3 errors listed in Table 1 are computed 
based on this assumption. As for the Mbh estimates, rever- 
beration mapping and stellar velocity dispersion methods 
give reliable estimates within factors of a few (Woo & Urry, 
2002). However, it is difficult to determine the actual error 
of all the individual Mbh estimates listed in Table 1. For 
that reason they will not be considered hereafter. 

Fig. 1 shows the alms vs Mbh plot. Clearly, there is 
a strong correlation between the two variables in the sense 
that that a rms decreases with increasing Mbh- Computa- 
tion of Kendal's r yields —0.69, which implies that the 
anti-correlation between al ma and Mbh is highly significant 
(probability > 99.5%). 

Let us suppose that the power spectrum, normalised to 
the mean squared, follows the form: P{y) = A{v /vhfb) -1 
Hz -1 , for v < fhft, and the form: P{v) = A(v/v h f b )~ 2 
Hz -1 , at higher frequencies (A is the PSD value at v = 
i>hfb)- Let us also suppose that the PSD "amplitude", i.e. 
PSD amp =Ax Uhfb, is constant and does not depend on Mbh, 
while p h f b is inversely proportional to the black hole mass, 
i.e. v hfb = C/M 7 Hz, where M 7 =M B h/1O 7 M . Under these 
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Table 1. The excess variance (<rj: ms ), PSD break frequency 
( v hfb) and black hole mass (Mbh) estimates of the objects stud- 
ied in this work. The objects are ranked by their Mbh- 



Name 


Type 


2 

^rms 


Vhfb 


BHMass 






(xl0~ 2 ) 


(xlO" 6 Hz) 


(xlO 7 M ) 


Akn 120 


BLS1 


3.7 ±0.7 


_ 


18.6(1) 


PG0804 


BLS1 


2.1 ±0.4 




17.6(1) 


3C120 


BLS1 


3.9 ±0.8 




13.5(2) 


NGC 5548 


BLS1 


5.4 ± 1.1 


0-6to' 9 5 (l) 


10.9(1) 


Fairal 9 (F9) 


BLS1 


4.4 ±0.9 


o-4tg:£(i) 


8.2(1) 


NGC 5506 


NLS1 




51.2tf . 8 (2) 


8.8* 


NGC 3516 


BLS1 


8.7 ± 1.7 


2.0^(1) 


1.7(3) 


NGC 4151 


BLS1 


10.3 ±2.1 


i-3±i:g(i) 


1.4(1) 


NGC 3783 


BLS1 


5.3 ± 1.1 


4.oi;-°(i) 


1.0(1) 


Mrk766 


NLS1 




500+ 3 3 °°(4) 


0.35(2) 


MCG-6-30-15 


NLS1 


8.0 ± 1.6 


100lJ°°(3) 


0.1(4) 


NGC 4051 


NLS1 


20 ±4.0 


800^(5) 


0.05(5) 


NGC 4395 


BLS1 




320+f 6 °(6) 


0.0066(6) 


Ark564 


NLS1 




1700l 6 9 °°(7) 


< 0.8(7) 



Column 2 lists the object type (BLS1 or NLS1). Note that 
NGC 5506 was recently classified as NLS1 (Nagar et al., 2002). 
All CT^ ms estimates in Col. (3) arc taken from Markowitz & 
Edelson (2001), except from PG 0804+761 which is taken from 
Papadakis et al. (2003). The numbers in parenthesis in Col. (4) 
correspond to the references for Vhfb as follows: (1) Markowitz 
ct al. (2003), (2) Uttley et al. (2002), (3) Vaughan et al. (2003), 
(4) Vaughan & Fabian (2003), (5) McHardy et al. (2003), (6) 
Shih et al. (2003), (7) Papadakis ct al. (2002). Note that the 
Vhfb values of NGC 5548 and F9, together with their errors, 
are taken from Table 5 of Markowitz et al., while the respective 
estimate of NGC 3783 was taken from Section 4.3 of the same 
work. The numbers in parenthesis in Col. (5) correspond to the 
references for M BH as follows: (1) Kaspi ct al. (2000), (2) Woo & 
Urry (2002), (3) Onken et al. (2003), (4) Uttley et al. (2002), (5) 
Shemmer et al. (2003), (6) Filippenko & Ho (2003), (7) Collier et 
al. (2001). The BH mass for NGC 5006 is estimated as explained 
in the text. Note that the M B h estimate of MCG 6-30-15 is 
rather uncertain (see discussion in Section 6.3 of Uttley et al., 
2002.) 




(Solar Mass) 



Figure 1. RXTE excess variance plotted as a function of Mbh- 
The solid line shows the best-fitting model in the case of a uni- 
versal PSD with u^fb oc Mbh -1 - The long-dashed line shows the 
"rms vs Mbh relation in the case when Vhfb scales with Mbh 
with a normalisation 20 times larger than case of the model shown 
with the solid line. The short-dashed and dot-dashed lines show a 
model with two PSD break frequencies, Vhfb and vif b , with vif b 
being 10 — 30 time smaller than Vhf b , respectively (see text for 
details). 



shown in Fig 1 (solid line) . The best fitting model parameter 
values are PSD amp = 0.017 ± 0.006, and C = 1.7±g;| x 10~ 6 
Hz (errors correspond to the 68% confidence for two in- 
teresting parameters, i.e. A\ 2 = X 2 + 2.3). Although the 
model does not provide a statistically accepted fit to the 
data (x = 14.8/6 dof), it describes rather well the overall 
trend of decreasing o 2 ms with increasing Mbh- The of ms 
values of the two NLSls (NGC 4051 and MCG 6-30-15), 
NGC 3783, and PG 0804±761 show the largest deviations 
from the best fitting model. 



assumptions, one can find a unique relation between a 2 ms 
and M7 using the fact that a 2 ms = J P(v)dv, where v\ s 
is the lowest frequency sampled (= 1/300 days = 3.858 x 
10~ 8 Hz, for the objects considered in this work). Note that, 
strictly speaking, the upper limit of the o 2 ma integral should 
be closer to 2AT obs , where AT obs ~ 1 ksec is the length 
of the individual monitoring exposures, since variations on 
shorter time scales are smoothed out. In any case though, 
use of either limits does not affect the integral's value almost 
at all. 

The integral's estimation depends on whether Uhfb is 
higher or lower than vi s . In the first case, one can show 
that: 

a 2 ms = PSD amp x [ln(C/M 7 ) + 18.07], (1) 

while in the second case, <J 2 ms is given by the relation, 

a 2 rma = PSD amp x 2.592 x 10 7 x C/M 7 . (2) 

The data shown in Fig. 1 were fitted with the model 
given by equations (1) and (2), excluding the points which 
correspond to the two NLSls. The best fitting model is also 



2.1 The possibility of a second frequency break in 
the PSDs 

If there is a second PSD flattening to zero slope at a fre- 
quency vifb<Vhfb (like the PSD of Galactic black hole can- 
didates, GBHs, when at low/hard state), then a 2 ma is given 
by the following relation, 

= PSD amp x [2 - vis I vi S b + ln(v hfb ) - ln(vifb)], (3) 

in the case when vifb> vi s . In the opposite case, u 2 ms 
is not affected by the presence of vif b and equations (1), 
(2) still hold. The dashed and dot-dashed lines in Fig. 1 
show a plot of equation (3) using the best fitting PSD amp 
and C values from the previous section, and assuming that 
Vhfb/vifb= 10 — 30, respectively. This ratio is similar to the 
ratio of the respective break frequencies in the PSD of Cyg 
X-l at low/hard state (Belloni & Hasinger, 1990). Equa- 
tion (3) provides an acceptable fit to all data (except from 
NGC 4051): x 2 = 10-7/7 dof and X 2 = 10-2/7 in the case 
when v h f b /vif b — 10 — 30, respectively. 

As Fig. 1 shows, the MCG 6-30-15 and NGC 3783 o 2 rms 
measurements are now consistent with the model defined 
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by equation (3). In other words, their af ms measurements 
suggest that there are two breaks in their PSD, both of 
which are higher than vi„. Interestingly, the observed PSD 
of NGC 3783 does show two breaks, with Vhfb/vifb^ 20 
(Markowitz et al., 2003). However, only one break (which 
corresponds to the "—2 to —1" slope flattening) has been 
reported in the case of MCG 6-30-15. Uttley et al. (2002) and 
Vaughan, Fabian, & Nandra (2003), report v h f b ~ 5 x 10~ 5 
Hz, and ~ 1 x 10~ 4 Hz, respectively, with the two values 
being consistent within the errors. If indeed Vhfb ~ 5 x 10~ 5 — 
10~ 4 Hz, the present results suggest that the low frequency 
break should be located at a frequency 10 — 30 times lower, 
i.e. at around ~ 1 — 5 x 10~ 6 Hz. Such a break should be 
detected by Uttley et al. (2002), as their PSD extends to 
frequencies lower than ~ 10~ 7 Hz. A combined long-term 
RXTE and short-term XMM-Newton PSD could resolve this 
issue. 



2.2 The case of NGC 4051 

NGC 4051 is the only source which deviates significantly 
from all model lines plotted in Fig. 1. This NLS1 is one 
of the most variable radio-quiet AGN. The discrepancy be- 
tween the NGC 4051 a^ m3 estimate and the model defined 
by equation (3) (dashed and dot-dashed lines in Fig. 1) sug- 
gests that the PSD of this object should have just one fre- 
quency break at frequencies higher than 1/300 days -1 . This 
is consistent with the results of McHardy et al. (2003) . Us- 
ing RXTE and XMM-Newton data, McHardy et al. detect a 
frequency break at ~ 8 x 10~ 4 Hz (where the PSD changes 
slope from —2 to —1), but they do not observe any further 
flattening at lower frequencies down to 10 -8 Hz. However, 
the NGC 4051 a rmg measurement does not agree with the 
best-fitting model defined by equations (1) and (2) (solid line 
in Fig. 1). This implies that either the PSD amp of NGC 4051 
is larger than in other AGN, or u h f b is higher than expected 
for its Mbh- According to the results of McHardy et al. 
(2003), PSD amp ~ 0.01, which is smaller than the best fit- 
ting PSD amp value of 0.017. Therefore, the a^ ma measure- 
ment of NGC 4051 can be consistent with equations (1) and 
(2) only if C is increased. The long-dashed line in Fig. 1 
shows a plot of the model defined by equations (1) and (2), 
using PSD am p= 0.017, and C = 20 times larger than the 
best-fitting value reported in Section 2. The agreement now 
between the model and the NGC 4051 measurement is good. 



2.3 The scaling of Vhfb with Mbh 

One of the results from the model fitting of the (tLs'Mbh 
relation is that, under the assumption of a universal PSD 
shape in AGN, then v h{b ~ 1.7 x 10~ 6 /M 7 Hz. It is inter- 
esting to see how this relation compares with the results 
from the power spectral analysis of recent observations in 
AGN. Table 1 lists the Mbh and Vhfb estimates for all AGN 
where a "—2 to —1" frequency break has been detected re- 
cently, irrespective of whether the same object has a long 
term £rj mj measurement or not. The Mbh of NGC 5506 was 
estimated using the stellar velocity dispersion measurement 
of 180 km/sec (Oliva et al, 1999), together with the Mbh - 
stellar velocity dispersion relation of Tremaine et al. (2002). 
Fig. 2 shows a plot of Uh fb as a function of Mbh • The 
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Figure 2. PSD Break frequency plotted as a function of Mbh- 
The solid line shows a plot of the v b f = C/M7 relation, with C 
being equal to the value from the best model fitting of the <J^ ms - 
Mbh relation. The dash line shows a plot of the same relation 
with C being 20 times larger. 

solid line is not a plot of the best fitting model to the data. 
Instead, it shows a plot of the v b f = C/M7 relation, with 
C = 1.7 x 10~ 6 Hz (i.e. the value estimated from the best 
model fitting to the ofms vs Mbh data in Section 2). The 
agreement between this line and the BLS1 data is very good. 
However, the NLSls fhfb estimates lie consistently above 
the solid line. The dashed line in Fig. 2 shows a plot of 
the v bf = C/M 7 with C = 20 x 1.7 x 10~ 6 Hz (i.e. the 
C value which can explain the high oYms measurement of 
NGC 4051, as discussed in Section 2.2). As Fig. 2 shows, the 
agreement between this line and the Vhfb measurements of 
NLSls is good. Therefore, both the large value of the <rj? ms 
measurement of NGC 4051 and the location of the NLS1 
data in Fig. 2 suggest that, for the same Mbh, Vhfb is ~ 20 
times higher in NLSls than in BLSls. 

2.4 A universal PSD amp value in AGN? 

Another result from the study the cr^ ms -MBH relation is that 
PSD omp should be ~ 0.017 in all AGN. It is not straightfor- 
ward to compare this result with recent results from power 
spectral analysis of X-ray light curves, as in most cases the 
PSD normalisation is not reported. In the cases where the 
relevant information is listed, the estimated PSD amp values 
are indeed close to 0.017. For example, using the results of 
Markowitz et al. (2003), one can find that: PSD amp =0.021, 
0.016, 0.012, 0.016, and 0.017 in the case of F9, NGC 5548, 
NGC 3783, NGC 3516, and NGC 4151, respectively. In the 
case of Cyg X-l, the Belloni & Hasinger (1990) data suggest 
that PSD amp = 0.02±0.002. This is consistent, within the er- 
rors, with the value of 0.017. Interestingly, the existing data 
for NLSls are systematically lower than 0.017. For exam- 
ple, Papadakis et al. (2002) and Vaughan & Fabian (2003) 
find PSD amp ~ 0.005 and ~ 0.01 in the case of Ark 564 and 
Mrk 766, respectively. Furthermore, McHrady et al. (2003) 
find PSD amp ~ 0.01 for NGC 4051. However, this value was 
estimated using a PSD model different than the "broken 
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Figure 3. The excess variance measurements plotted as a func- 
tion of the X— ray (2 — 10 kcV) luminosity. The dot-dashed line 
shows a plot of the expected c^ ms — Lx relation using the results 
from the best model fitting of the <rJ: ms -MBH relation (shown 
with the dot-dashed line in Fig. 1), and assuming that Lx ocMgjj- 
The solid and dashed lines show plots of models defined to ac- 
count for the high cr^ ms measurement of NGC 4051 (see text for 
details). 

power- law model" used in the other cases. In fact, as Fig- 
ure 17 in McHardy et al.(2003) shows clearly, the PSD amp 
of NGC 4051 is probably as large as the PSD amp of other 
BLSls. 



2.5 The accretion rate of AGN 

It is interesting to investigate if the of ms vs 2 — 10 keV 
X-ray luminosity (Lx) correlation of Markowitz & Edelson 
(2001) is in fact a byproduct of the <J^ ms vs Mbh correlation 
presented in this work. This is possible only if there exists a 
relation between Mbh and Lx , which is the same for all the 
objects. This implies that the accretion rate, when consid- 
ered as a fraction of the Eddington limit, should be similar 
in all objects. 

Fig. 3 shows the cr$ ms -Lx relation for the objects with 
o^ms measurements listed in Table 1, using the Lx mea- 
surements of Markowitz et al. (2003) for F9, NGC 5548, 
NGC 3783, NGC 3516, and NGC 4151, Markowitz & Edel- 
son (2001) for 3C120, Aknl20, and MCG -6-30-15, McHardy 
et al. (2003) for NGC 4051, and the measurement of Pa- 
padakis et al. (2003) for PG 0804+761. The dot-dashed line 
in Fig. 3 shows a plot of a model, which is based on: (1) the 
best fitting "2 PSD breaks with isifb/fhfb=30" model to the 
ff?ms-MBH relation (the one shown with the dot-dashed line 
in Fig. 1), and (2) the relation L x = (BH Mass /M s ) x 10 35 ' 8 
ergs/sec. As Fig. 3 shows clearly, this line fits the data very 
well (x 2 = 8.8/9 dof), except from NGC 4051. Using the 
mean bolometric conversion factor of L bo i ~ 27 x Lx of 
Padovani & Rafanelli (1988), the adopted relation between 
L x and Mbh implies that L bo i = 1.7 x lO 37 (BH M a SS /M ) 
ergs/sec. Since Lsdd = 1-25 x lO 38 (BHMass/M0) ergs/sec, 
the main conclusion is that the <j^ ma -Lx correlation could 
be explained in terms of the uLs-Mbh relation if all the 



objects in the present sample radiate at ~ 13.5% of their 
Eddington luminosity. 

The fact that NGC 4051 is not consistent with the dot- 
dashed line plotted in Fig. 3 implies that either its accretion 
rate is significantly larger than ~ 10% of the Eddington 
limit, or its o% ma measurement is significantly higher than 
that expected for its luminosity. The solid line in Fig. 3 
shows the expected al ma -Lx relation, using: (1) the best 
fitting "1 PSD break with fhfb ot Mbh" 1 " model to the 
u ms -MBH relation (shown with the solid line in Fig. 1) 
and (2) assuming an accretion rate three times higher than 
0.1 of the Eddington limit. The agreement between the 
model and NGC 4051 is reasonably good. The dashed line 
in Fig. 3, shows the expected <j rms -Lx relation, using: (1) 
the best fitting "1 PSD break with fhfb oc Mbh 1 and C 
20 times larger" model to the crJ ms -MBH relation (dashed 
line in Fig. 1) and (2) assuming an accretion rate ~ 0.1 
of the Eddington limit. This model also agrees well with 
NGC 4051. Therefore, the "unusually" large variability am- 
plitude of NGC 4051 for its luminosity does not necessarily 
imply an exceptionally high accretion rate; a PSD break at 
a frequency higher than what is expected for its Mbh can 
also account for the discrepancy of NGC 4051, with respect 
to the other objects in the sample. 



3 DISCUSSION AND CONCLUSIONS 

The main results of this work are the following: 

(i) There is a significant correlation between the long term 
excess variance of AGN and their Mbh, in the sense that 
u rms decreases with increasing Mbh. 

(ii) The <r rms -MBH relation is consistent with the hypothesis 
of a universal PSD shape in AGN with two frequency breaks, 
fhfb and vifb- The PSD changes its slope from —2 to —1 at 
fhfb, and from —1 to zero at vifb- 

(iii) The high frequency break, fhfb, decreases with increas- 
ing Mbh as v h fb= 1.7 x 10~ 6 Mf 1 Hz. The low frequency 
break, vifb, is ~ 10 to 30 times smaller than Vhfb- 

(iv) The PSD peak amplitude in the v x P(y) space is ~ 
0.017 in all objects. 

(v) The correlation between <r 2 ms and Lx is a byproduct of 
the o-rm S -MBH correlation, with AGN radiating at ~ 10 — 
15% of the Eddington luminosity. 

(vi) The a^ms measurement of NGC 4051, a NLS1 object, 
is larger than what is expected for its Mbh and Lx- This 
result can be explained if its PSD has only one break (at 
frequencies higher than ~4x 10~ 8 Hz) which is located at 
a frequency ~ 20 times higher than the frequency expected 
in the case of a BLS1 with the same Mbh- An intrinsically 
higher PSD amplitude (e.g. see McHardy et al. 2003) may 
also contribute to the larger crj? ms of this source. 

These results are based on ofm S measurements of long 
light curves in the 2—10 keV band. Since the power spec- 
tra of individual AGN appear to be energy dependent (e.g. 
Papadakis & Lawrence 1995, Nandra & Papadakis 2001, 
Vaughan et al. 2003), if there is indeed a "universal" PSD 
shape, this should also be energy dependent. For example, 
a slope of —2 should characterize the AGN power spectra 
above the high frequency break in the 2 — 10 keV band only. 

The results listed above are entirely consistent with the 
recent results from the power spectral analysis of RXTE, 
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XMM-Newton, and ASCA light curves. For example, in all 
cases where a high frequency PSD slope break has been de- 
tected, the slope of the PSD below and above this break is 
roughly equal to —1 and —2 (see e.g. Uttley et al., 2002; 
Markowitz et al., 2003). Furthermore, as discussed in Sec- 
tion 2.3, the v/i/b-MaH relation determined from the cr^ ms - 
Mbh correlation is entirely consistent with the existing data. 
In fact, the relation of T break (= l/i/ h/6 )=M B H/(10 6 ' 5 M ) 
days, of Markowitz et al. (2003), is in agreement with the 
results of the present work. Finally, the observed PSD amp 
values are also consistent with the value of 0.017, which is 
derived from the study of the uLs-Mbh relation. 

One of the most interesting results presented in this 
work is that the 2-10 keV band PSD of BLSls, at least, 
has a "universal" shape: there two break frequencies, the 
slope between them is ~ —1, the slope above the high fre- 
quency break ~ — 2, the high frequency break is mass depen- 
dent, but the PSD amp is not. Instead it has a constant value 
of ~ 0.02. Since this result is based on the measurement 
of a^ms an d Mbh of ten objects only (most of which are 
among the most frequently observed AGN in X-rays) it is 
then premature to accept that it holds for the Seyfert galax- 
ies as a class. Even if it does, it is not expected that either 
the high frequency slope slope for example or PSD amp will 
be exactly equal to —2 and 0.02 for each individual AGN. 
Instead, these values should be considered as "average" or 
"typical" values around which the respective values of all 
AGN will be distributed, in the same way that the energy 
spectral slopes of AGN are distributed around the "canon- 
ical" value of r ~ 1.9. Confirmation of this result can be 
achieved in two ways. First, by using the results from the 
PSD analysis of many AGN. This is the direct and most 
powerful method, as it will allow us to actually determine 
the distribution of the PSD slope and PSD amp values around 
the "canonical" values. However, its rather unlikely that the 
number of AGN with well studied PSDs will increase sig- 
nificantly in the near future. What the present works shows 
clearly is that there is a second way, which involves the al ma 
estimation of many objects with known Mbh- Although the 
CTrms °f an individual object cannot constrain its PSD shape 
in any way, the study of the it™ s -Mbh relation can be useful 
in order to investigate whether there is indeed a "universal" 
PSD shape in AGN. Furthermore, contrary to the former 
possibility, the use of archival data of many objects for the 
detailed study of the of ms -MBH relation should be possible 
in the near future as discussed in the Introduction. 

If indeed the PSD of BLSls AGN has a "universal" 
shape (in the sense that was described above) then this re- 
sult could constrain physical models proposed to explain the 
physical process responsible for the X-ray variations. For ex- 
ample, models should be able to predict the i^/6-Mbh re- 
lation, and explain the "universal" PSD amp value of ~ 0.02 
in AGN (perhaps in GBHs as well). As an example of how 
this can be achieved, I present below a simple "exercise", 
based on the assumption that the X-ray variations of AGN 
at frequencies between vi fb and Vh fb are associated with ac- 
cretion disc instabilities which cause variations either to the 
energy release in the innermost regions of the disc or to the 
soft photons input to the X-ray emitting corona. 

As a first step, one should investigate the relevance of 
physical processes to the origin of X-ray variability by com- 
paring the various physical time-scales in accretion disc to 



the PSD break time-scales. Let us consider a BLS1 with 
Mbh= 10 7 Mq. According to the i/|,/5-Mbh relation, the 
break time-scale for this object should be Tbreak = 6.8t^g 
days. Let us also consider the case of a standard a— disc 
(Shakura & Sunyaev, 1973) around the central object, and 
let us assume that Tbreak is associated with one of the 
characteristic time-scales of the disc, i.e. the orbital (T orb ), 
thermal (T t h), sound-crossing (T sc ), and viscous time-scale 
(T V isc). Using the relations given by Treves, Maraschi, & 
Abramowicz (1988), one can estimate these time-scales at 
every radius from the central object. The results show 
that Tbreak ~ T or b at R ~ 70 — SORs (where Rs is the 
Schwarzschild radius). However, this seems rather unlikely 
as most of the energy in the accretion disc is dissipated 
at smaller radii. Tbreak could also correspond to T sc at 
R ~ 3 — 5Rsch (assuming that R ~ 100//, where H is the 
disc's scale height). However, if T br eak for a NLS1 with the 
same Mbh is ~ 20 times faster, and T brea k corresponds to 
the same physical time-scale in both type of objects, then 
obviously this cannot be the sound-crossing time-scale. Fur- 
thermore, T visc is far too slow to account for Tbreak- As- 
suming that a < 0.1 (where a is the viscosity parameter), 
and /?<//, then T vi3C > 870 days even at R = 3Rs- Fi- 
nally, Tbreak could correspond to T t h at R ~ 30 — AQRseh, 
for a = 0.1. In the case of a NLS1 with the same Mbh, 
Tbreak could also correspond to T t h, but at a smaller radius 
(~4-6-R s ). 

Viscosity variations which happen at large radii, and 
develop on the local T visc , can affect the energy release of the 
accretion disc in its innermost part, and result in a PSD with 
slope —1 between ui/t and Vhfb (Lyubarskii, 1997). However, 
if the PSD breaks are associated with the thermal time- 
scales, then perhaps as the accretion rate increases at various 
radii and the disc becomes unstable, thermal instabilities 
set in. The disc temperature will increase, increasing the 
local flux, hence the soft photons input to the X-ray corona. 
These soft photon "flares" should then develop on the local 
thermal time-scale. Assuming an exponential increase of soft 
photons with time, the superposition of these flares with the 
different time-scales can lead to a PSD with a flat shape up 
to vifb= l/T th (R out ), a -1 slope up to vifb= l/T th (Ri n ), 
and —2 slope at higher frequencies (Lehto 1989), where Rout 
and Rin are the outer and inner radius of the disc region 
which is affected by the thermal instabilities. If R ou t = 200 — 
300/? s , and R m = 30 - 40/?s, then the resulting PSD will 
be in agreement with the present results. 

The fact that the P{v) x v value between v\ f b and v h f b is 
similar in all objects implies that the PSD integral between 
vijb and Vhfb (i.e. the excess variance associated with the 
components of those frequencies) remains roughly constant, 
irrespective of Mbh- In the context of shot noise models 
where flares occur randomly with a constant rate of say A 
flares per unit time, and have the same shape and amplitude, 
then of ms oc 1/A. Considering the case of soft photons flares, 
which develop on the local thermal time-scale, as Mbh in- 
creases, and T t h increases accordingly, it is natural to expect 
that A should decrease, as it should take longer time for the 
instability to develop and the accretion disc to "relax" before 
the next instability starts to build up. In this case we should 
expect PSDam P to increase with increasing Mbh. However, 
at the same time, as Mbh increases, the flares should become 
"smoother" (because Tth also increases). This fact should re- 
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suit in a decrease in 07 ms . Perhaps then, because an increase 
of the Mbh causes a decrease of A and, at the same time, 
"smoother" flares, and because these two effects affect the 
observed variability in the opposite way, ofms mav turn out 
to be roughly constant in all objects. However, further work 
is needed in order to examine if such a model can also pre- 
dict correctly the PSD amp value of ~ 0.02, as suggested by 
this work. 

The differences in the variability properties between 
NLSls and BLSls can be explained by the fact that Vhfb 
is higher in a NLS1 than in a BLS1 with the same Mbh. 
The same explanation has also been proposed recently by 
McHardy et al. (2003). In the framework of the picture pre- 
sented above, which involves thermal instabilities at a cer- 
tain region of the accretion disc, the different v h f b values 
should correspond to differences in the location of this region 
in the two type of objects. In the innermost region of BLSls 
(i.e. at radii smaller 30Rs) the thermal instabilities could 
be stabilized by a physical mechanism (e.g. the formation of 
a hot wind). The same mechanism may not be sufficient to 
stabilize the innermost region of the NLSls, perhaps because 
the accretion rate is larger in these objects. Whether NLSls 
show a second PSD break at vif b ~ 10 — 30 times smaller 
than Vhfb is not certain yet. The results of this work show 
that the long term a^ ms measurement of NGC 4051 is not 
consistent with the presence of a second break, in agreement 
with the results of McHardy et al.(2003). As these authors 
point out, NGC 4051 is the first Seyfert galaxy which shows 
a PSD similar to the power spectrum of GBHs at high/soft 
state. However, the cr^ ms measurement of MCG -6-30-15 is 
consistent with the presence of a second break. Ark 564 also 
shows two breaks in its PSD (Pounds et al., 2001; Papadakis 
et al., 2002). Power spectrum analysis of more NLSls or the 
study of the drms-MBH relation of a large number of NLSls 
is necessary to clarify this issue. 
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